Adiabatic loading of one-dimensional SU(iV) alkaline earth fermions in optical lattices 
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Ultracold fermionic alkaline earth atoms confined in optical lattices realize Hubbard models with 
internal SU(A^) symmetries, where A'^ can be as large as ten. Such systems are expected to harbor 
exotic magnetic physics at temperatures below the superexchange energy scale. Employing quan- 
tum Monte Carlo simulations to access the low-temperature regime of one-dimensional chains, we 
show that after adiabatically loading a weakly interacting gas into the strongly interacting regime 
of an optical lattice, the final temperature decreases with increasing TV. Furthermore, we esti- 
mate the temperature scale required to probe correlations associated with low-temperature SU(A'^) 
magnetism. Our findings are encouraging for the exploration of exotic large- A'^ magnetic states in 
ongoing experiments. 

PACS numbers: 67.85.-d,03.75.Ss,37.10. Jk 



Ultracold fermionic alkaline earth atoms confined in 
optical lattices realize an important, tunable generaliza- 
tion of the Hubbard model, widely used to model strongly 
correlated electrons [H-Ql- Within this generalization, 
the conventional spin-1/2 SU(2) symmetry is enhanced 
to an SU(A^) symmetry. = 2/ -|- 1 is determined by 
the nuclear spin. /, which varies in alkaline earths from 
/ = 1/2 to / = 9/2 depending on the atomic species. 
The SU(A^) symmetry arises because the electronic de- 
grees of freedom have neither spin nor orbital angular 
momentum (due to the closed-shell structure) and thus 
decouple from the nuclear spin. In the electronic ground 
state the symmetry has been theoretically predicted to 
hold to an accuracy of 10~^ [ij and experiments have 
constrained deviations to be less than 5 x 10~^ De- 
spite the large /, quantum fluctuations remain important 
due to the enhanced symmetry, giving rise to magnetic 
frustration and exotic ground states, such as valence- 
bond solids M- 111, exotic spin ordcrings [lij or (chiral) 
spin liquids |13l . Il4l | in addition to the plethora of poten- 
tial phases for the conventional N = 2 Hubbard model, 
such as antiferromagnets, d-wave superconductors, and 
nematic states. The possibility of mimicking such excit- 
ing many-body p hysics, as well as potential applications 
to atomic clocks [l5Hl8|, measurements of fundamental 
constants [l^ , and quantum information processinj 
has stimulated substantial experimental progress [; 

Although cold atom experiments routinely reach 
nanokelvin temperatures, it is an ongoing effort to 
achieve temperatures and entropies sufficiently low to see 
superexchange driven magnetic many-body physics 0, 
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33j . Hazzard et al. showed via a high temperature se- 
ries expansion (HTSE) that for final temperatures T >t, 
with t the timneling rate in the lattice, the temperatures 



reached after adiabatic loading from experimentally real- 
istic initial conditions decreases with increasing N 34 1. 



Closely related is the finding that the entropy of Mott 
states increases with N even faster than the initial weakly 
interacting gas entropy. The question of the behavior 
below t, and especially below the magnetic exchange en- 
ergy scale J ^ t'^/U where U is the on-site interaction, 
has remained open even though this is one relevant for 
exploring exotic SU(A^) magnetism. Two issues are par- 
ticularly relevant in this regime: how does A^ affect (i) 
the temperature reached by adiabatic loading and (ii) the 
physical properties, such as correlation functions? 

Here, we address both questions in one dimensional 
systems using quantum Monte Carlo calculations. We 
show that also for T < t'^/U the temperatures reached 
by adiabatic loading decrease with increasing A^. This 
decrease occurs even relative to the temperature scales 
of interesting physics, for example the onset of Luttinger 
liquid behavior, magnetic correlations, or ground state- 
like correlations. 

The Hamiltonian of the SU(A^) Hubbard model de- 
scribing alkaline earth atoms in optical lattices is [l[ 



n = -t 
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(1) 



where /J, ^ {fa,i) arc creation (annihilation) operators for 

fermions of flavor a at site i, and nf = f^^fa.i- At filling 
1/A^ (density n = 1) local moment (Heisenberg) quantum 
magnetism arises in the strong interaction limit, where 
second order processes lead to a flavor exchange interac- 
tion scale J ~ /U . We employ quantum Monte-Carlo 
(QMC) to obtain the low temperature thermodynamic 
properties of one-dimensional lattice alkaline earth atoms 
described by Eq. ^ in the relevant regime T <^t ^U, 
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FIG. 1: (Color online) Left: Final vs. initial temperature 
for the adiabatic loading scheme for A'^ = 2, 3, 4, 6 and 10 
(from left to right) inside a harmonic trap with trap frequency 
uj — 2n X 90Hz and for U/t — 8. One can alternatively trans- 
late the initial temperature into units of Fermi temperature 
Tf, using Ti/Tp = {2.8,3.2,3.5,4.1,4.8} x IQ-^iTi/uj) for 
N = {2,3,4,6,10}. The inset shows the temperature de- 
crease relative to the N — 2 Hubbard model, rjv' = Tf{N — 
2)/Tf{N = iV'), for (from bottom to top) TV' = 3 (red), 
iV' = 4 (green), iV' = 6 (blue) and TV' = 10 (magenta) vs. 
Ti/u}. Right: Distance measure Ag of the magnetic correla- 
tions relative to the ground state values for a homogeneous 
system with U /t = 8 ai n = 1. 



much lower than the temperature scales accessible within 
the HTSE. 

These thermodynamic properties allow us to calculate 
the temperatures that experiments can achieve using the 
standard experimental adiabatic lattice loading protocol: 
the systems are prepared by starting with a weakly in- 



teracting gas and slowly turning on the lattice [32, 
Ideally, this is done sufficiently slowly to maintain adi- 
abaticity, minimizing the increase in entropy. The final 
temperature in the adiabatic limit in this sense provides 
a lower bound to the experimentally achievable temper- 
atures, and in practice the adiabatic bound is frequently 
an accurate approximation (sil . [ssl . [35| . 

In the following, we consider loading a three dimen- 
sional {d = 3) trapped gas (no lattice) into a two di- 
mensional array of independent one dimensional lattices. 
Fig. [T] (left panel) shows our central result, the final tem- 
perature Tf of the harmonically trapped lattice system 
with fixed trapping frequency w as a function of the 
initial temperature Ti before applying the lattice [36| . 
We use realistic experimental particle number and trap 
frequencies, similar to Ref. [3J]: M = 1.5 x 10^ and 
u; = 27r X 90Hz. For loading ^'^^yb in a 266 nm lat- 
tice, our choice of U corresponds to the lattice depth 
required to obtain U/t = 8. We observe from Fig. [T] 
that the final temperature achieved for Ti/w < 9 signif- 
icantly decreases with increasing TV. This includes the 
lowest temperatures, an order of magnitude lower than 



where the HTSE is applicable, and is encouraging for ex- 
perimentally achieving SU(TV) quantum magnetism. For 
instance, at the relatively warm temperature Ti/uj = 4 
(some experiments are already at even lower tempera- 
tures) the final temperature is already decreased by fac- 
tors of 1.7 (TV = 3), 1.85 (TV = 4), 3.25 (TV = 6) and 5.6 
(TV = 10) compared to the conventional TV = 2 case. We 
note that even this is a pessimistic estimate since these 
ratios compare different TV fixing Ti/uj. Experimentally, 
Ti will likely decrease with TV, as discussed in Ref. |34. 
The reason why Tf decreases with TV, for fixed initial T^, 
is the scaling of the initial and final states' entropy with 
TV. At very low temperature T < /U , the Mott in- 
sulator possesses TV — 1 gapless channels, and the metal 
. The increasing number of gapless 



possesses TV [37 



excitations implies that the entropy growth at fixed tem- 
perature is faster than the initial state's N^^^. Thus, 
the final temperatures will decrease with increasing TV. 
This is even more favorable than at high temperatures 
T > t, where the Mott insulator accommodates an en- 
tropy S ^ log TV as can be seen in Fig. [T] However, 
even for T > t, for TV < 20 the logarithmic term grows 
fast enough to compensate for the N^^^ in the initial 
state, leading to colder final states with increasing TV. 

The enhanced cooling effect has in deed been observed 
(at elevated temperatures) for SU(6) fermions in recent 
experiments on a ^^'^Yb gas in a three-dimensional opti- 
cal lattice (39} . In particular, their findings are in good 
agreement with the HTSE results (valid for their exper- 
iments), confirming the validity of the adiabaticity as- 
sumption for their experiments. 

We now describe the procedure used to obtain Fig. [H 
Assuming adiabaticity, the initial temperature Ti and 
particle number Mi in the absence of the lattice uniquely 
determine the final temperature Tf through the conserva- 
tion of entropy and particle number. Tf can be calculated 
in terms of the experimentally measurable Ti and A/i once 
one knows how the particle number and entropy of the 
system in a lattice depend on temperature T and the 
chemical potential /i, J\ff{T,fi) and Sf{T,^), as follows: 
We calculate Tf and ^/ by solving particle number and 
entropy conservation equations, Nf{Tf,^f) = Mi(Ti,^i) 
and Sf(Tf,^f) = Si(Ti,^i), for Tf and /z/ in terms of 
the initial temperature Ti and chemical potential /i^ of 
the weakly interacting gas. In practice, since experi- 
mentalists measure the initial particle number A/i rather 
than ^i, we solve for /x^ in terms of Tj and Mi using 
Ni{^i,Ti) ~ Ni- The conservation equations become 

^^f{Tf,^Jif)^K, Sf{Tf,fif)^S,{T,,fi.,iM.,T,)). (2) 

The functions J\fi{T,^) and Si{T,iJ,) are those of a gas 
in a harmonic trapping potential V{r) = mw^r^/2. For 
relevant experimental initial conditions, the temperature 
is low compared to the Fermi temperature, interactions 
are weak, and the number of particles is large, so the 
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We con- 



semiclassical approximation to the non-interacting de- 
generate Fermi gas is accurate, giving [s^] Ni{n,T) 
iN/d\){fi/ujy emdS.ifi.T) 
sidcr dimension d = 3. 

We obtain the functions Mf{T, fj.) and Sf{T, fi) in two 
steps. First, to compute the total particle and entropy 
in the trap we apply the local density approximation 
(LDA) [40|: we obtain the properties at position r from 
those of a homogeneous system at a chemical potential 
^(r) = ^ — V{r). Due to the large particle number (about 
30 particles in the central tube at T = 0), this approxima- 
tion will be accurate. Second, to obtain the homogeneous 
system properties used in the LDA, we use sign-problem- 
free QMC simulations 4l| within the stochastic series 



expansion (SSE) framework [42| - |44| . We first calculate 
the density and entropy for finite systems up to L = 100 
sites, with open boundaries within the grand-canonical 
ensemble for various values of the chemical potential /i. 
We find finite size effects to be negligible on these lattice 
sizes for the quantities of interest below. The entropy 
is obtained by a standard thermodynamic integration of 
the energy E. In fact, the HTSE agrees well with the 
QMC data down to T/t w 2, as illustrated in Fig. d 
and we thus used the HTSE results at T/t = 10 as the 
high-temperature end point for the thermodynamic in- 
tegration. Increasingly dense temperature grids are re- 
quired to perform this integration down to lower temper- 
atures, and we were able to perform this procedure down 
to T/t = 0.1. We find that our data indeed connects to 
the T limit, where 5" scales linear in T [s!]- We fi- 
nally compared our QMC results to ground state proper- 
ties (ssj . obtained using the density matrix renormaliza- 
tion group (DMRG) |45l448{ . In particular, the extrap- 
olated ground state energies agree with DMRG results 
within the statistical error bars. Finally, we calculate 
tables of A// and Sf for a dense grid of /x and T, con- 
struct interpolating functions of this data supplemented, 
for very negative /x, with the virial expansion that ac- 
curately describes the gas in this regime. The adiabatic 
loading equations Eqs. ^ are then solved numerically to 
obtain the results in Fig. [1] 

Having determined achievable temperatures, we next 
turn to correlation functions that are indicative for low- 
energy SU(A^) magnetism in these systems. The relevant 
spin-spin correlation function is related to the equal and 
unequal flavor density-density correlation functions [s^ 



Si\^-j\) 



-E 



N{N - 1) 



This can be measured via, for example, Bragg spec- 
troscopy [i^ . To compare finite-temperature and ground 
state correlations (obtained from DMRG), we define an 
appropriate distance measure that accounts also for the 




FIG. 2: (Color online) Left: Entropy S per site for iV = 3, 
IJ /t = 8 at fixed chemical potentials from pt/t = —2 (bottom) 
to fi/t — 8 (top) obtained from QMC for a homogeneous 
system. We note, that the system forms a Mott insulating 
ground state with density n = 1 for ^/t ~ 1.5 to 5. Dashed 
hues down to T/t = 1 denote HTSE results from Ref. HI. 
For clarity, an offset was added to the entropies, indicated 
by the arrows. Right top: Temperature dependence of the 
energy per site E relative to the extrapolated ground state 
energy Eo, from QMC. Right middle: Temperature depen- 
dence of the specific heat Cv [13] • The arrow indicates the 
low-T peak in Cv in the Mott regime. Right bottom: Tem- 
perature dependence of the kinetic energy Ekin, from QMC. 
Statistical errors are below the line width. 



overall magnitude of the ground state correlations [38|, 

This quantity is sensitive to both long range and short 
range magnetic correlations (sol . [sH . The right panel of 
Fig. [T] shows A5 for n = 1 and [//i = 8 as a function of 
Ti/w, using Tf/t determined from the lattice loading [H^ . 
This shows that increasing N brings the system not only 
to lower values of T, but also into a region where ground- 
state-like correlations are more developed. 

A hallmark of quantum magnetism in A^-componcnt 
Luttingcr liquids (LL) are 2kF = 27rn/7V oscillations in 
the spin correlation functions, with n the density [ssllssj. 
Figure Fig. [3] shows the spin structure factor <S(fc), the 
Fourier transformation of iS(r), for A'^ = 3 and 4 at 
U /t = 8 and fixed filling of 1/A^ (n = 1). At low tempera- 
tures T/t < 0.1, the 2k F peak appears, as expected from 
LL theory. This is associated with a maximum in the 
specific heat Cy [HJl shown in Fig. [51 marking the onset 
of SU(Af) Heisenberg physics [HI]. At high temperatures 
T/t > 1, S{k) instead becomes featureless. However, in 
an intermediate regime 0.5 ^ T/t < 1, a broad peak at 
A: = TT is observed, which shifts towards 2/cf as the tem- 
perature is lowered towards T t^ /U: at these temper- 
atures, nearest neighbor correlations emerge. While this 
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FIG. 3: (Color online) Spin structure factor S{k) for a homo- 
geneous system with density n = 1 for TV = 3 (left) and 
4 (right), U/t = 8 as a function of the final temperature 
T/t (here, L = 48). The characteristic 2fcF peaks below 
T jt ~ 0.15 signal the onset of magnetic correlations. 



regime is not representative for the ground state behav- 
ior, the short-ranged magnetic correlations determine the 
main features of the structure factors, and experiments 
reaching this temperature regime should be able to in- 
vestigate magnetic behavior. In Fig. [21 the entropy and 
kinetic energy similarly display the onset of the mag- 
netic exchange region as a kink and large decrease, re- 
spectively. This happens around T jt k, 0.5. For N = 2^ 
we have checked that our data are consistent with former 
studies of the SU(2) Hubbard model 

For intermediate temperature /U '^T t, one can 
ask whether the system realizes a SU(-/V) generalization 



of the spin- incoherent Luttinger liquid (siLL) j57l [58 1 
where the charge degress of freedom show LL-like be- 
haviour with a simultaneous absence of significant spin 
correlations. Here, we find no clear signatures of siLL in 
the density structure factor near k = im (not shown). 
We leave a more in-depth study for future work, where 
it will likely be beneficial to examine the momentum dis- 
tribution n(fc), as was done for the SU(2) case [ssf . 

Summary. — We studied the one-dimensional SU(iV) 
Fermi-Hubbard model at finite temperatures employing 
numerically exact quantum Monte Carlo simulations. We 
calculated the density and entropy as functions of chem- 
ical potential and temperature and used this to deter- 
mine final temperatures of the lattice system after adi- 
abatically loading an optical lattice from a degenerate 
gas. We found substantial decreases of the final temper- 
ature with increasing N, even down to low temperatures 
T < t'^/U <^ t. Together with our results for the temper- 
ature dependence of correlation functions, we envisage 
that it should be possible to explore features of SU(7V) 
quantum magnetism in ongoing experiments with alka- 
line earth atoms. 
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Note added. — Recently, Messio and Mila explored 
the A'^-dependence of the entropy and correlations in 
the Heisenberg limit of the SU(A^) systems considered 
here (Hoj : their results are in accord with our finding for 
the Hubbard model, obtained at U/t = 8. 
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